Reversible transitions between noradrenergic and mesenchymal tumor identities define cell plasticity in neuroblastoma.

Noradrenergic and mesenchymal identities have been characterized in neuroblastoma cell lines according to their epigenetic landscapes and core regulatory circuitries. However, their relationship and relative contribution in patient tumors remain poorly defined. We now document spontaneous and reversible plasticity between the two identities, associated with epigenetic reprogramming, in several neuroblastoma models. Interestingly, xenografts with cells from each identity eventually harbor a noradrenergic phenotype suggesting that the microenvironment provides a powerful pressure towards this phenotype. Accordingly, such a noradrenergic cell identity is systematically observed in single-cell RNA-seq of 18 tumor biopsies and 15 PDX models. Yet, a subpopulation of these noradrenergic tumor cells presents with mesenchymal features that are shared with plasticity models, indicating that the plasticity described in these models has relevance in neuroblastoma patients. This work therefore emphasizes that intrinsic plasticity properties of neuroblastoma cells are dependent upon external cues of the environment to drive cell identity.

Single-cell transcriptomics reveals shared immunosuppressive landscapes of mouse and human neuroblastoma.

BACKGROUND: High-risk neuroblastoma is a pediatric cancer with still a dismal prognosis, despite multimodal and intensive therapies. Tumor microenvironment represents a key component of the tumor ecosystem the complexity of which has to be accurately understood to define selective targeting opportunities, including immune-based therapies. METHODS: We combined various approaches including single-cell transcriptomics to dissect the tumor microenvironment of both a transgenic mouse neuroblastoma model and a cohort of 10 biopsies from neuroblastoma patients, either at diagnosis or at relapse. Features of related cells were validated by multicolor flow cytometry and functional assays. RESULTS: We show that the immune microenvironment of MYCN-driven mouse neuroblastoma is characterized by a low content of T cells, several phenotypes of macrophages and a population of cells expressing signatures of myeloid-derived suppressor cells (MDSCs) that are molecularly distinct from the various macrophage subsets. We document two cancer-associated fibroblasts (CAFs) subsets, one of which corresponding to CAF-S1, known to have immunosuppressive functions. Our data unravel a complex content in myeloid cells in patient tumors and further document a striking correspondence of the microenvironment populations between both mouse and human tumors. We show that mouse intratumor T cells exhibit increased expression of inhibitory receptors at the protein level. Consistently, T cells from patients are characterized by features of exhaustion, expressing inhibitory receptors and showing low expression of effector cytokines. We further functionally demonstrate that MDSCs isolated from mouse neuroblastoma have immunosuppressive properties, impairing the proliferation of T lymphocytes. CONCLUSIONS: Our study demonstrates that neuroblastoma tumors have an immunocompromised microenvironment characterized by dysfunctional T cells and accumulation of immunosuppressive cells. Our work provides a new and precious data resource to better understand the neuroblastoma ecosystem and suggest novel therapeutic strategies, targeting both tumor cells and components of the microenvironment.

BET and CDK Inhibition Reveal Differences in the Proliferation Control of Sympathetic Ganglion Neuroblasts and Adrenal Chromaffin Cells.

Neuroblastoma arising from the adrenal differ from ganglionic neuroblastoma both genetically and clinically, with adrenal tumors being associated with a more severe prognosis. The different tumor properties may be linked to specific tumor founder cells in adrenal and sympathetic ganglia. To address this question, we first set up cultures of mouse sympathetic neuroblasts and adrenal chromaffin cells. These cultures were then treated with various proliferation inhibitors to identify lineage-specific responses. We show that neuroblast and chromaffin cell proliferation was affected by WNT, ALK, IGF1, and PRC2/EZH2 signaling inhibitors to a similar extent. However, differential effects were observed in response to bromodomain and extraterminal (BET) protein inhibitors (JQ1, GSK1324726A) and to the CDK-7 inhibitor THZ1, with BET inhibitors preferentially affecting chromaffin cells, and THZ1 preferentially affecting neuroblasts. The differential dependence of chromaffin cells and neuroblasts on BET and CDK signaling may indicate different mechanisms during tumor initiation in sympathetic ganglia and adrenal.

ALK mutation dynamics and clonal evolution in a neuroblastoma model exhibiting two ALK mutations.

The ALK gene is a major oncogene of neuroblastoma cases exhibiting ALK activating mutations. Here, we characterized two neuroblastoma cell lines established from a stage 4 patient at diagnosis either from the primary tumor (PT) or from the bone marrow (BM). Both cell lines exhibited similar genomic profiles. All cells in the BM-derived cell line exhibited an ALK F1174L mutation, whereas this mutation was present in only 5% of the cells in the earliest passages of the PT-derived cell line. The BM-derived cell line presented with a higher proliferation rate in vitro and injections in Nude mice resulted in tumor formation only for the BM-derived cell line. Next, we observed that the F1174L mutation frequency in the PT-derived cell line increased with successive passages. Further Whole Exome Sequencing revealed a second ALK mutation, L1196M, in this cell line. Digital droplet PCR documented that the allele fractions of both mutations changed upon passages, and that the F1174L mutation reached 50% in late passages, indicating clonal evolution. In vitro treatment of the PT-derived cell line exhibiting the F1174L and L1196M mutations with the alectinib inhibitor resulted in an enrichment of the L1196M mutation. Using xenografts, we documented a better efficacy of alectinib compared to crizotinib on tumor growth and an enrichment of the L1196M mutation at the end of both treatments. Finally, single-cell RNA-seq analysis was consistent with both mutations resulting in ALK activation. Altogether, this study provides novel insights into ALK mutation dynamics in a neuroblastoma model harbouring two ALK mutations.

Study of chromatin remodeling genes implicates SMARCA4 as a putative player in oncogenesis in neuroblastoma.

In neuroblastoma (NB), genetic alterations in chromatin remodeling (CRGs) and epigenetic modifier genes (EMGs) have been described. We sought to determine their frequency and clinical impact. Whole exome (WES)/whole genome sequencing (WGS) data and targeted sequencing (TSCA®) of exonic regions of 33 CRGs/EMGs were analyzed in tumor samples from 283 NB patients, with constitutional material available for 55 patients. The frequency of CRG/EMG variations in NB cases was then compared to the Genome Aggregation Database (gnomAD). The sequencing revealed SNVs/small InDels or focal CNAs of CRGs/EMGs in 20% (56/283) of all cases, occurring at a somatic level in 4 (7.2%), at a germline level in 12 (22%) cases, whereas for the remaining cases, only tumor material could be analyzed. The most frequently altered genes were ATRX (5%), SMARCA4 (2.5%), MLL3 (2.5%) and ARID1B (2.5%). Double events (SNVs/small InDels/CNAs associated with LOH) were observed in SMARCA4 (n = 3), ATRX (n = 1) and PBRM1 (n = 1). Among the 60 variations, 24 (8.4%) targeted domains of functional importance for chromatin remodeling or highly conserved domains but of unknown function. Variations in SMARCA4 and ATRX occurred more frequently in the NB as compared to the gnomAD control cohort (OR = 4.49, 95%CI: 1.63-9.97, p = 0.038; OR 3.44, 95%CI: 1.46-6.91, p = 0.043, respectively). Cases with CRG/EMG variations showed a poorer overall survival compared to cases without variations. Genetic variations of CRGs/EMGs with likely functional impact were observed in 8.4% (24/283) of NB. Our case-control approach suggests a role of SMARCA4 as a player of NB oncogenesis.

Activated ALK signals through the ERK-ETV5-RET pathway to drive neuroblastoma oncogenesis.

Activating mutations of the ALK receptor occur in a subset of neuroblastoma tumors. We previously demonstrated that Alk mutations cooperate with MYCN overexpression to induce neuroblastoma in mice and identified Ret as being strongly upregulated in MYCN/Alkmut tumors. By a genetic approach in vivo, we now document an oncogenic cooperation between activated Ret and MYCN overexpression in neuroblastoma formation. We show that MYCN/RetM919T tumors exhibit histological features and expression profiles close to MYCN/Alkmut tumors. We show that RET transcript levels decrease precedes RET protein levels decrease upon ALK inhibition in neuroblastoma cell lines. Etv5 was identified as a candidate transcription factor regulating Ret expression from murine MYCN/Alkmut tumor transcriptomic data. We demonstrate that ETV5 is regulated both at the protein and mRNA levels upon ALK activation or inhibition in neuroblastoma cell lines and that this regulation precedes RET modulation. We document that ALK activation induces ETV5 protein upregulation through stabilization in a MEK/ERK-dependent manner. We show that RNAi-mediated inhibition of ETV5 decreases RET expression. Reporter assays indicate that ETV5 is able to drive RET gene transcription. ChIP-seq analysis confirmed ETV5 binding on the RET promoter and identified an enhancer upstream of the promoter. Finally, we demonstrate that combining RET and ALK inhibitors reduces tumor growth more efficiently than each single agent in MYCN and AlkF1178L-driven murine neuroblastoma. Altogether, these results define the ERK-ETV5-RET pathway as a critical axis driving neuroblastoma oncogenesis downstream of activated ALK.

Heterogeneity of neuroblastoma cell identity defined by transcriptional circuitries.

Neuroblastoma is a tumor of the peripheral sympathetic nervous system, derived from multipotent neural crest cells (NCCs). To define core regulatory circuitries (CRCs) controlling the gene expression program of neuroblastoma, we established and analyzed the neuroblastoma super-enhancer landscape. We discovered three types of identity in neuroblastoma cell lines: a sympathetic noradrenergic identity, defined by a CRC module including the PHOX2B, HAND2 and GATA3 transcription factors (TFs); an NCC-like identity, driven by a CRC module containing AP-1 TFs; and a mixed type, further deconvoluted at the single-cell level. Treatment of the mixed type with chemotherapeutic agents resulted in enrichment of NCC-like cells. The noradrenergic module was validated by ChIP-seq. Functional studies demonstrated dependency of neuroblastoma with noradrenergic identity on PHOX2B, evocative of lineage addiction. Most neuroblastoma primary tumors express TFs from the noradrenergic and NCC-like modules. Our data demonstrate a previously unknown aspect of tumor heterogeneity relevant for neuroblastoma treatment strategies.

HMCan-diff: a method to detect changes in histone modifications in cells with different genetic characteristics.

Comparing histone modification profiles between cancer and normal states, or across different tumor samples, can provide insights into understanding cancer initiation, progression and response to therapy. ChIP-seq histone modification data of cancer samples are distorted by copy number variation innate to any cancer cell. We present HMCan-diff, the first method designed to analyze ChIP-seq data to detect changes in histone modifications between two cancer samples of different genetic backgrounds, or between a cancer sample and a normal control. HMCan-diff explicitly corrects for copy number bias, and for other biases in the ChIP-seq data, which significantly improves prediction accuracy compared to methods that do not consider such corrections. On in silico simulated ChIP-seq data generated using genomes with differences in copy number profiles, HMCan-diff shows a much better performance compared to other methods that have no correction for copy number bias. Additionally, we benchmarked HMCan-diff on four experimental datasets, characterizing two histone marks in two different scenarios. We correlated changes in histone modifications between a cancer and a normal control sample with changes in gene expression. On all experimental datasets, HMCan-diff demonstrated better performance compared to the other methods.

Activated Alk triggers prolonged neurogenesis and Ret upregulation providing a therapeutic target in ALK-mutated neuroblastoma.

Activating mutations of the ALK (Anaplastic lymphoma Kinase) gene have been identified in sporadic and familial cases of neuroblastoma, a cancer of early childhood arising from the sympathetic nervous system (SNS). To decipher ALK function in neuroblastoma predisposition and oncogenesis, we have characterized knock-in (KI) mice bearing the two most frequent mutations observed in neuroblastoma patients. A dramatic enlargement of sympathetic ganglia is observed in AlkF1178L mice from embryonic to adult stages associated with an increased proliferation of sympathetic neuroblasts from E14.5 to birth. In a MYCN transgenic context, the F1178L mutation displays a higher oncogenic potential than the R1279Q mutation as evident from a shorter latency of tumor onset. We show that tumors expressing the R1279Q mutation are sensitive to ALK inhibition upon crizotinib treatment. Furthermore, our data provide evidence that activated ALK triggers RET upregulation in mouse sympathetic ganglia at birth as well as in murine and human neuroblastoma. Using vandetanib, we show that RET inhibition strongly impairs tumor growth in vivo in both MYCN/KI AlkR1279Q and MYCN/KI AlkF1178L mice. Altogether, our findings demonstrate the critical role of activated ALK in SNS development and pathogenesis and identify RET as a therapeutic target in ALK mutated neuroblastoma.

Breakpoint features of genomic rearrangements in neuroblastoma with unbalanced translocations and chromothripsis.

Neuroblastoma is a pediatric cancer of the peripheral nervous system in which structural chromosome aberrations are emblematic of aggressive tumors. In this study, we performed an in-depth analysis of somatic rearrangements in two neuroblastoma cell lines and two primary tumors using paired-end sequencing of mate-pair libraries and RNA-seq. The cell lines presented with typical genetic alterations of neuroblastoma and the two tumors belong to the group of neuroblastoma exhibiting a profile of chromothripsis. Inter and intra-chromosomal rearrangements were identified in the four samples, allowing in particular characterization of unbalanced translocations at high resolution. Using complementary experiments, we further characterized 51 rearrangements at the base pair resolution that revealed 59 DNA junctions. In a subset of cases, complex rearrangements were observed with templated insertion of fragments of nearby sequences. Although we did not identify known particular motifs in the local environment of the breakpoints, we documented frequent microhomologies at the junctions in both chromothripsis and non-chromothripsis associated breakpoints. RNA-seq experiments confirmed expression of several predicted chimeric genes and genes with disrupted exon structure including ALK, NBAS, FHIT, PTPRD and ODZ4. Our study therefore indicates that both non-homologous end joining-mediated repair and replicative processes may account for genomic rearrangements in neuroblastoma. RNA-seq analysis allows the identification of the subset of abnormal transcripts expressed from genomic rearrangements that may be involved in neuroblastoma oncogenesis.

Characterization of rearrangements involving the ALK gene reveals a novel truncated form associated with tumor aggressiveness in neuroblastoma.

Activating mutations of the ALK gene have been identified in sporadic and familial cases of neuroblastoma (NB), a cancer of the peripheral nervous system, and are thought to be the primary mechanism of oncogenic activation of this receptor in this pediatric neoplasm. To address the possibility that ALK activation may occur through genomic rearrangements as detected in other cancers, we first took advantage of high-resolution array-comparative genomic hybridization to search for ALK rearrangements in NB samples. Using complementary experiments by capture/paired-end sequencing and FISH experiments, various types of rearrangements were fully characterized, including partial gains or amplifications, in several NB cell lines and primary tumors. In the CLB-Bar cell line, we described a genomic rearrangement associated with an amplification of the ALK locus, leading to the expression of a 170 kDa protein lacking part of the extracellular domain encoded by exons 4 to 11, named ALK(Δ4-11). Analysis of genomic DNA from the tumor at diagnosis and relapse revealed that the ALK gene was amplified at diagnosis but that the rearranged ALK allele was observed at the relapse stage only, suggesting that it may be implicated in tumor aggressiveness. Consistently, oncogenic and tumorigenic properties of the ALK(Δ4-11) variant were shown after stable expression in NIH3T3 cells. Moreover, we documented an increased constitutive kinase activity of this variant, as well as an impaired maturation and retention into intracellular compartments. These results indicate that genomic rearrangements constitute an alternative mechanism to ALK point mutations resulting in receptor activation.

The CDKN2A/CDKN2B/CDK4/CCND1 pathway is pivotal in well-differentiated and dedifferentiated liposarcoma oncogenesis: an analysis of 104 tumors.

The MDM2 and CDK4 genes are the main targets of chromosome 12 amplification in well-differentiated and dedifferentiated liposarcomas. Nevertheless, around 10% of these tumors do not amplify CDK4. To find substitutive alterations of CDK4 amplification, we analyzed a large series of liposarcomas by array-CGH, real-time genomic PCR, gene expression array, and real-time RT-PCR. We demonstrate that an alteration in the CDKN2A/CDKN2B/CDK4/CCND1 pathway is present in almost all cases without CDK4 amplification, thereby confirming the pivotal role of this pathway in liposarcoma oncogenesis. Moreover, we show that cell cycle and differentiation are driven by a subtle and complex balance between members of this pathway. Finally, we demonstrate that in tumors without amplification/overexpression of CDK4, the chromosome 1q21-1q23 region is a preferential partner of chromosome 12 amplicon, suggesting that the mechanism of amplification is slightly different in this group of tumors.

MDM2 and CDK4 immunohistochemistry is a valuable tool in the differential diagnosis of low-grade osteosarcomas and other primary fibro-osseous lesions of the bone.

Low-grade osteosarcoma is a rare malignancy that may be subdivided into two main subgroups on the basis of location in relation to the bone cortex, that is, parosteal osteosarcoma and low-grade central osteosarcoma. Their histological appearance is quite similar and characterized by spindle cell stroma with low-to-moderate cellularity and well-differentiated anastomosing bone trabeculae. Low-grade osteosarcomas have a simple genetic profile with supernumerary ring chromosomes comprising amplification of chromosome 12q13-15, including the cyclin-dependent kinase 4 (CDK4) and murine double-minute type 2 (MDM2) gene region. Low-grade osteosarcoma can be confused with fibrous and fibro-osseous lesions such as fibromatosis and fibrous dysplasia on radiological and histological findings. We investigated MDM2-CDK4 immunohistochemical expression in a series of 72 low-grade osteosarcomas and 107 fibrous or fibro-osseous lesions of the bone or paraosseous soft tissue. The MDM2-CDK4 amplification status of low-grade osteosarcoma was also evaluated by comparative genomic hybridization array in 18 cases, and the MDM2 amplification status was evaluated by fluorescence in situ hybridization or quantitative real-time polymerase chain reaction in 31 cases of benign fibrous and fibro-osseous lesions. MDM2-CDK4 immunostaining and MDM2 amplification by fluorescence in situ hybridization or quantitative real-time polymerase chain reaction were investigated in a control group of 23 cases of primary high-grade bone sarcoma, including 20 conventional high-grade osteosarcomas, two pleomorphic spindle cell sarcomas/malignant fibrous histiocytomas and one leiomyosarcoma. The results showed that MDM2 and/or CDK4 immunoreactivity was present in 89% of low-grade osteosarcoma specimens. All benign fibrous and fibro-osseous lesions and the tumors of the control group were negative for MDM2 and CDK4. These results were consistent with the MDM2 and CDK4 amplification results. In conclusion, immunohistochemical expression of MDM2 and CDK4 is specific and provides sensitive markers for the diagnosis of low-grade osteosarcomas, helping to differentiate them from benign fibrous and fibro-osseous lesions, particularly in cases with atypical radio-clinical presentation and/or limited biopsy samples.